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Abstract 

We address some recent erroneous claim that Hq observations are difficult to accommodate with 
LTB cosmological models, showing how to construct solutions in agreement with an arbitrary value 
of Hq by re- writing the exact solution in terms of dimensionless parameters and functions. This 
approach can be applied to fully exploit LTB solutions in designing models alternative to dark energy 
without making any restrictive or implicit assumption about the inhomogeneity profile. The same 
solution can also be used to study structure formation in the regime in which perturbation theory is 
not enough and an exact solution of the Einstein's equation is required, or to estimate the effects of 
a local inhomogeneities on the apparent equation of state of dark energy. 
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I. INTRODUCTION 



Inhomogeneous solutions of the Einstein's field equations have received same attention as 
alternatives cosmological model to dark energy. The main idea is that the space inhomogeneity 
is able to mimic in red-shift space the effects of homogeneous dark energy, without the need 
of modifying gravity or introducing a model for dark energy. Some interesting theoretical 
investigation about LTB solutions in a cosmological context can be found for example in 



24| . Another interesting application of LTB solutions is to consider them as realistic model of 



the non-linear effects of large scale inhomogeneities in presence of dark energy 
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, and recently 



25|-|32|. There 



37] it has been 



have been several attempts to fit cosmological data 
claimed that accurate measurements of Hq are sufficient to rule out best fit void models. In 
this paper we will show how such claims are erroneous, and it is easy to construct general LTB 
models which can accommodate any observed value of Hq value, arguing that the large residual 
functional freedom should be fully exploited before making conclusive statements about their 
compatibility with observations. Our approach differ from previous ones because Hq is not 
a quantity derived from other parameters of the model, but rather the fundamental natural 
physical scale in terms of which all other dimensionful quantities are expressed. The paper is 
organized as following. Using a generalized conformal time variable we derive in a unified and 
compact way the solution with and without cosmological constant , showing how the latter one 
reduces to the first one in the appropriate limit. We then introduce dimensionless functions 
defining the LTB solution, expressing the analytical solution in a form in which the observed 
Hq is a manifestly free length scale parameter. The analytical solutions we derive could have 
different applications; they could for example be be used to fit observational cosmological data 
or for an exact study of the spherical collapse in structure formation. 



II. EXACT SOLUTION FOR LTB MODELS 



The LTB solution can be written as 
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4l| as 



+ R^dQ'^ 



1) 



l + 2^(r 

where i? is a function of the time coordinate t and the radial coordinate r, E{r) is an arbitrary 
function of r, and R r = drR(t, r). 
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The Einstein equations with dust and a cosmological constant give 

(R\' _ 2E{r) , 2M(r) , A 

p^'^'^ = wr;/ 

with M(r) being an arbitrary function of r, i? = dtR{t, r) and c = SvrG = 1 is assumed 
throughout the paper. Since Eq. ([2]) contains partial derivatives respect to time only, its 
general solution can be obtained from the FLRW equivalent solution by making every constant 
in the latter one an arbitrary function of r. 

An analytical solution can be found by introducing a new coordinate r] = rj{t, r) and a 
variable a by 

( drj\ r _ 1 



dt R a 



(4) 



and new functions by 



Then Eq. ([2]) becomes 



_ 6M(r) _ 2E{r) 



where a is now regarded as a function of t] and r, a = a{r],r). It should be noted that the 
coordinate r], which is a generalization of the conformal time in a homogeneous FLRW universe, 
has been only implicitly defined by Eq. (jl]). The actual relation between t and r] can be obtained 
by integration once a{r], r) is known: 

t{T],r) = / a{x,r)dx + tb{r) , (7) 
Jo 



which can be computed analytically, and involve elliptic integrals of the third kind 

The function ^^(r) plays the role of constant of integration, and is an arbitrary function of r, 
sometime called bang function, since by construction at time t = th{r) we have a{tiy{r),r) = 0, 
and correspond to the fact that the big bang initial singularity can happen at different times 
at different positions from the center in a LTB space. 



A. Case with no cosmological constant 



When A = the eq.(l6]) can be easily solved, and the solution of physical interest satisfying 
the the big bang initial condition a(0, r) = is : 
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a{r],r) 



Po{r) 
6k{r) 

and after integration, using eq.(I 
t{ri,r) 



1 — cos ( \J k{r) rj 



po{r) sin^ {j\/k(r)v 
3k{r) 



(8) 



Po[r) 
6k{r) 



V 



k{r) 



sin I J k{r) rj 



+ hir) . 



(9) 



It should be noted that this solution is valid for any sign of the the function fc(r), and is 
more convenient than the equivalent two branches form used by some other authors [? ]. 



B. Case with cosmological constant 

The general analytical solution for a FLRW model with dust and cosmological constant was 



obtained by Edwards [42| in terms of elliptic functions. By an appropriate choice of variables 
and coordinates, we may extend it to the LTB case thanks to the spherical symmetry of both 
LTB and FLRW models, and to the fact that dust follows geodesies without being affected 
by adjacent regions. Inspired by the construction of the solution for the FLRW case and by 
dimensional analysis we can make the ansatz 

a(T,r) = ^ , (10) 

^ ' 6K{r) + -f(P{Ta) ^ ' 

where we are using the following dimensionless quantities: 

k{r) = (aoHoYKir) , (11) 

r] = T{aoHo)-\ (12) 

Po(r) = 3nl,ir)alH', , (13) 

A = S^aH^ , (14) 

a{r],r) = a{T{aoHoy\r) = a{T,r) , (15) 
da{r],r) da{T,r) rr \ 

— — — = aoHo — — — = aoHoa (T, r) , (16) 

and we are denoting with a tilde the function (p for consistency with the notation for the 
scale factor, i.e. we use a tilde for quantities which depend on the dimensionless generalized 
conformal time variable T. After substituting the above ansatz in the Einstein's equation we 
get a differential equation equation of the form: 

A4>'{Tf + B4){Tf + C4){Tf + D^{T) + E = Q, (17) 



where the coefficients {^4, B, C, D} depend on the parameters {a, (5, 7, 5}. 
If we now impose the conditions 

I - -4 (18) 

C = (19) 
^ = ^72 (20) 

\= 9z (21) 

the Einstein's equation f lT7|l reduces to the canonical form of the Weierstrass differential equa- 
tion 

/ ~ \ ^ 

(^j =40^-^20-^3, (22) 

whose solution is an elhptic function conventionally expressed using the notation (j){T, g2, g2)- 
The above conditions give : 

^-12^' 

(25) 

= 432^^ • ^26) 

Since we have only two conditions for the four unknowns {a, /3, 7, S} the ansatz we made allows 
to express the solution of the Einstein's equation in terms of the Weierstrass elliptic functions 
for an arbitrary choice of two of them. We choose : 

{a = l;P = 3aon%{r) ; 7 = 12 ; 5 = 1} , (27) 

which gives the solution 

~a(T r) = ^^o^lii^) .28) 

ir(r) + 120(T,^2(r),<73(r))' 

92{r) = ^, (29) 



\-{2Kirf-27Qj,{QUr)r) , (30) 



or after multiplying every term by (aoi^o)^ and using the original dimensionful quantities 

r],k{r),poir) 
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k{r) + 12(f){r],g2{r),g3{r)) 
(j){ri,r) = (j){r]{aoHo),r){aoHof = (j){T,r){aoHof . 



(32) 



C. Relation between the two solutions 

In the previous subsections we have derived the analytical solution of the Einstein's equa- 
tion for a LTB model with and without cosmological constant. Since we are using the same 
coordinate system to derive the two solutions we expect that eq. (!27|) should reduce to eq.([8]) 
when A = 0. In fact the discriminant of the Weierstrass function which solves eq.(l6]) 



is zero when the cosmological constant vanishes. 

It can be shown that when the discriminant is zero the Weierstrass elliptic function can be 
expressed in terms of trigonometric functions : 



If we now substitute the above equation in eq.( l27|) . after some manipulation we get exactly 
eq.([8]), which can be considered a consistency check for the solutions derived so far. In the rest 
of the paper we choose the so called FLRW gauge, i.e. the coordinate system in which po{r) is 
constant. 

III. LTB SOLUTIONS AND Hq OBSERVATIONS 

The exact solutions derived above are not defined in terms of directly observables quantities, 
since we have not introduced explicitly any system of units. If we are interested in cosmological 
applications the natural scale to choose is Hq^, so it is convenient to re-write equation (JH]) in 
the familiar FLRW-like form: 



A = gl- 21 gl = ±^nj,(nUr)f {AK{rf - 27 , n^{nUr)f) 



(33) 




(34) 




(35) 



(36) 



where we have define 



H[t,r) 


a 
a 


(37) 


TT 

-no 


TJ ( + (W 


(38j 








Kir) 




(40) 


ilK[t,r) 


= -^('■)(ir) • 


(41) 




= ii°MW(^)', 


(42) 


Po(r) 


= 3Q^(r)aX^ 


(43) 


A 




(44) 



In this form the observed Hq appears exphcitly as one of the parameters defining the LTB 
model, and all the other quantities are dimensionless. 

The exact solution can then be conveniently re-written in terms of dimensionless functions 
and coordinates as 



a{T,r) 



aol^^(r) sin^ iTJir(r 



K{r) 



t{T,r) = 



° 2K{r) 
rj = T{aoHo)-'. 



T - 



4W) 



sin J Kir) T 



+ U{r) , 



After defining the Hubble rate in terms of the generalized conformal time variable r] as 



a{t,r) a{r]^rY 



air, r) 



(45) 

(46) 

(47) 
(48) 

(49) 
(50) 



if we want the exact solution to be in agreement with Hq observations we can impose the two 
following conditions 



a(77o,0) = ao, 
^''^''(^0,0) = Ho, 



(51) 

(52) 



where Oq is, as expected, an arbitrary parameter, and r/o is the value of the generalized conformal 
time coordinate r] corresponding to the central observer today, i.e. satisfying 



t{r]Q, 0) = to 
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(53) 



We will discuss the consequences of the above conditions in the two following sections according 
to the presence or absence of the cosmological constant term in the Einstein's equations. As 
mentioned at the end other previous section, we will will use the freedom in the choice of the 
radial coordinate to fix the coordinates in which po = const., and consequently Q%j{r) = const. 



IV. VANISHING COSMOLOGICAL CONSTANT CASE 



In the vanishing cosmological constant case we get 

KoifS + 1) " 



ao 



^0 , 



where 



K{0) 



To = VoiaoHo 



fo 



tan 



(54) 
(55) 

(56) 
(57) 
(58) 



and we have kept ao, Hq on both side of the equations just to show their relation with eqs. fl51|55p . 
but they obviously cancel. The above equations can be solved for fo, to get 



fo 



1 + Ko 



or equivalently, after inverting eq.( l58|l 



"A/ 



To 



l + Ko, 

2arctan(/o) 2 arctan (y/Ko) 



(59) 
(60) 

(61) 
(62) 



VKo VKo 

The meaning of the above relations is that for an arbitrary value of Ko there is always a LTB 
solution in agreement with i^o observations given by 



a(r, r) 



ao(Ko + l)sin2 {\T^Kiy) 



K{r) 



tiT, r) = 
To 



il + Ko 



2K{r 
2 arctan (^/Ko 



T 



K{r) 



sin JK{r) T 



^0 



(63) 
(64) 

(65) 



As expected ao does not appear in observable quantities such as the cosmic time t{7],r). This 
solution will be by construction in agreement with an arbitrary value of Hq, while the residual 
functional freedom in K{r) could be exploited to fit other observables such as the luminosity 
distance. We will report about this in a separate work, in particular in relation to previous failed 
attempts to fit both the luminosity distance and the cosmic microwave radiation (CMB). What 
we have found is in agreement with a previous investigation based on a local Taylor expansion 



43|, and provides a useful method to design LTB models without cosmological constant to fit 



observational data. 



V. NOT VANISHING COSMOLOGICAL CONSTANT CASE 

The case in which the cosmological constant is not vanishing can be interesting because it 
provides a simple way to model spherically symmetric large scale inhomogeneities around a 
central observer, which, as discussed in are expected to be the dominant ones due to other 
strong evidences of isotropy such as the CMB radiation. After re-expressing the metric using 
the same dimensionless quantity introduced above we get: 



ao , (66) 



00 - ^0, (67) 
where 

0o = 0(To,O), (68) 

- (i) . («^) 

/ {T=To,r=0} 

Ko = K{0). (70) 

Our final goal is to determine Qm^Tq from the above equation, but apparently we have an extra 
unknown, i.e. 0q. This is not really independent from (pQ, since the Weierstrass equation allows 
to express algebraically one in terms of the other. This is a general property which can be used 
whenever we have to simplify or manipulate analytical expressions involving the solution we 
have derived in terms of elliptic functions: 



oyo ^ AV MJ ^^^3^ ^j^-^ 



216 12 16 
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Finally the solution of the system of equations. (151II55P gives: 



nl = Ko-Qa + I, (72) 
4>o = ^(2Ko-3fiA + 3), (73) 

Another possible way to interpret the solutions we have derived is to consider the Einstein's 
equation at the center which gives 

i = nl{o) + nK + nA, (74) 

where we are using the notation 

Qk = -K{0) , (75) 

since locally at the center the LTB solution is equivalent of a FLRW solution with the above 
value of Qk- This is in fact equivalent to the equation fl55|) . since the solution of the Einstein's 
equation substituted into the definition of H^^^ in eq. (l49l) gives the Einstein's equation itself. 
According to this interpretation the Einstein's equation at the center gives the condition to 
determine in agreement with the Hq observations. Using the equations derived above, for 
a given value of f^A, Kq we can now determine and 4>q which will give a model in agreement 
with an arbitrary observed Hq. We can finally solve eq.f l69|) to get Tq or rjQ. While in the 
vanishing cosmological constant case this can be done analytically, in this case we need to solve 
numerically the equation to get Tq. 

VI. POSSIBLE APPLICATIONS 

One of the main problems that LTB models without dark energy seem to face in fitting 



observational data is the value of Hq 35|, l36|]. The approach normally adopted previously 



consisted in defining a void model in terms of a limited set of parameters and then indirectly 
derive the value of Hq. In other words Hq was not considered as a fundamental scale in terms 
of which to determine the geometry of the LTB solution as we are proposing here. In our 
approach instead Hq is the primary ingredient in defining the LTB model, and this should help 
to improve the fitting of these model to observations. In order to fit CMB data a LTB model 
should have the same distance to the last scattering surface of the best fit AC DM model, and 
have the same value of H{zls) at the last scattering surface. It is also known from a numerical 
inversion method [? ], that it is possible to construct a LTB model able to exactly mimic both 
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the luminosity distance Dl{z) and the H{z) of a AC DM model, but this requires the use of 
both k{r) and ti,{r). This implies that allowing for a inhomogcncous big bang we can expect to 
fit both luminosity distance and the CMB, and in fact previous claims of the impossibility to 
fit CMB, Hq and Dl[z) were based on void models with ^^(r) = 0. Since using the solution in 
the form we have derived gives a LTB model in agreement with any value of Hq, the problem 
will probably be in getting the value of H{zls) in agreement with CMB observations, but a 
conclusive answer to this question can be given only by fitting the data with models which 
have more freedom in defining the inhomogeneity profile, while the void models analyzed so 
far e only depend on few parameters. In this regard a differential inversion method would be 
more appropriate since it would allow to fully exploit the functional degree of freedom of LTB 
solutions, his application is left to a separate work. The important point is that contrary to 
previous claims to rule out void models using Hq observations, the correct way to design a LTB 
model is to introduce Hq from the beginning and express the other parameters and functions 
such as K{r) in terms of it, rather than deducing it. The indirect method used so far in fitting 
cosmological data can in fact artificially limit the genuine freedom of LTB solutions according 
to the particular ansatz made for the inhomogeneity profile. 



VII. CONCLUSIONS 



We have shown how to construct a LTB solutions in agreement with any value of Hq by 
using it as the fundamental scale in terms of which to define all the other quantities determining 
the model, proving that any claims that Hq measurements are enough to rule out LTB models 
are wrong. This will allow to fit data in a more efficient way, since Hq does not need to be 
computed from the other parameters of the model, but on the contrary can be fixed directly 
without any additional computational effort. In this way it will be possible to fully explore the 
space of LTB solutions as viable cosmological model, going beyond the limited class of models 
examined so far. 

Another natural application is in relation to structure formation, where the derived solutions 
could be used to study exactly the spherical collapse in presence of a cosmological constant. 
In this context it could actually be more convenient to define the solution directly in terms of 
the density contrast, and we leave this to a future work. The exact solution could also be used 
in numerical simulations of structure formation in the regime in which the collapse of compact 
objects cannot be treated accurately with perturbation theory, and an exact solution of the 
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Einstein's equations is required. 
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